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function (ft but the total charge of the topological defects is still unchanged. 
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1. Introduction 



• i 

The world of topological defects is amazingly rich and have been the focus of much attention in 
many areas of contemporary physics ]l|, ^ [|. The importance of the role of defects in understanding 
a variety of problems in physics is clear 0, ^ || [/J. Recently, some physicists noticed^, |9| that 
the topological defects are closely related to the spontaneously broken of 0(m) symmetry group to 
0(m — 1) by m-component order parameter field <ft and pointed out that for m — 1, one has domain 
walls, m — 2, strings and m — 3, monopoles, for m = 4, there are textures. But for the lack of a 
powerful method, the topological properties are not very clear, the unified theory of describing the 
topological properties of all these defect objects is not established yet. 

In this paper, in the light of 0-mapping topological current theory |10[, a useful method which 
plays a important role in studying the topological invariant sflTT|, an< ^ ^ ne topological structures of 



physical systems |T3|, O, IT3J, we will investigate the topological quantization and the branch process 
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of arbitrary dimensional topological defects. We will show that the topological defects are generated 
from where = and are topological quantized under the condition J(-) 7^ 0. While at the zero 
points of field function <fi where the corresponding Jacobian determinant J(^) vanishes, the defect 
topological current bifurcates and the topological defects split or merge at such point, this means 
that the topological defects system is unstable at these points. 

This paper is organized as follows. In section 2, we investigate the topological quantization of 
these topological defect and point out that the topological charges of these defects are the Winding 
numbers which are determined by the Hopf indices and the Brouwer degrees of the 0-mapping. 
In section 3, we study the branch process of the defect topological current at the limit points, 
bifurcation points and higher degenerated points systematically by virtue of the 0-mapping theory 
and the implicit function theorem. 



2. Topological quantization of topological defects 

In our previous papers mentioned above, only the topological current of point-like particles was 
discussed. In this paper, in order to study the topological properties of arbitrary dimensional topo- 
logical defects, we will extend the concept to present an arbitrary dimensional generalized topological 
current. We consider the 0-mapping as a map between two manifolds, while the dimensions of the 
two manifolds are arbitrary. It is an important generalization of our previous work on topological 
current and is of great usefulness to theoretical physics and differential geometry. 

In n-dimensional Riemann manifold G with the metric tensor and local coordinates x M v = 
1, n), a m-component vector order parameter field <j)(x) can be looked upon as a mapping between 
the Riemann manifold G and a m-dimensional Euclidean space R m 

(f): G — > R m , <j) a = (j) a (x), o = l,.. 

The direction field of <fi(x) is generally determined by 

« a W = \\4>{x)\\ = v>(*)0 a (*) (i) 

with 

n a (x)n a (x) = 1. (2) 



It is obviously that n a (x) is a section of the sphere bundle S(G) [fLO] . If n a (x) is a smooth unit vector 
field without singularities or it has singularities somewhere but at the point <p(x) 7^ 0, from (H) we 
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have 

n a d^n a = 0, fi = l,...,n, (3) 

which can be looked upon as a system of n homogeneous linear equations of n a (a = 1, ...,m) with 
coefficient matrix [c^ra"] . The necessary and sufficient condition that (^) has non-trivial solution for 
n a {x) is rank [9 M n a ] < m, i.e. the Jacobian determinants 

D^(dn) = ^■■■^+^e ai ... am d tlh+1 n^ ■ ■ ■ d, n n a ™ (4) 

are equal to zero, where k = n — m. While, at the point = 0, the above consequences are not held. 
In short, we have the following relations 

D^{dn)[ = °' ^>°> (5) 
1/0, for = 0, 

which implies D' 11 '" flk (dn) behaves itself like a function 5(<f)). So we are focussed on the zeroes of 
<p a {x). 

Suppose that the vector field 0(x) possesses I isolated zeroes, according to the implicit function 
theorem|16[], when the zeroes are regular points of 0-mapping at which the rank of the Jacobian 



matrix [<9^0 a ] is m, the solutions of <fi = can be expressed parameterizedly by 

X» = Z?{u\--;U k ), t = l,...,l, (6) 

where the subscript i represents the i-th solution and the parameters u 1 (I = 1, ...,k) span a k- 
dimensional submanifold with the metric tensor gjj = g^v^j^p- which is called the i-th singular 
submanifold N{ in the Riemannian manifold G corresponding to the 0-mapping. For each singular 
manifold Ni, we can define a normal submanifold Mi in G which is spanned by the parameters v A 
with the metric tensor qab = Qhv^x^b (A, B = I,...,m), and the intersection point of Mj and 
iVj is denoted by Pi which can be expressed parameterizedly by v A = pf. In fact, in the words of 
differential topology, Mj is transversal to iVj at the point pi, i.e. 

T Pi (G)=T P XM l )+T Pi (N l ). 

By virtue of the implicit function theorem, it should be held true that, at the regular point Pi, the 
Jacobian matrices J(^) satisfies 

j A = Mi^n a o (7) 
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In the following, we will induce a rank-fc topological current through the integration of D^ 1 "'^ k {dn) 



in (||) on Mi. As is well known, the generalized Winding Number fl7H has been given by the Gauss 
map n : <9£; -> S™' 1 

W ' = ^(S-^-DI L, A • ■ • A (8) 

where 

A{S m ~ l ) - 



' T(m/2) 

is the area of (m — l)-dimensional unit sphere S 171 " 1 , n* denotes the pull back of map n and <9Ej the 
boundary of a neighborhood of p« on Mj with pi ^ <9Xj, Ej fl T,j = 0. The generalized Winding 
Numbers Wi can also be rewritten as 



W- = / e n ai dn a2 A • • • A dn am 

which means that, when the point or v A covers <9£j once, the unit vector n a will cover a region 
n[9£j] whose area is Wi times of A(S m ~~ l ), i.e. the unit vector n a will cover the unit sphere S 771 ^ 1 
Wi times. From the above equation, one can deduce that 

1 

AiS^Ym- 1)! J 9Ml 

= -T77, — ttt rr / ea,-a d„.„ n ai da h ^n a2 • • • d u n am dx» k+1 A • • • A dx» n 

A(S m ~ 1 )(m — 1)! J Mi 1 m ^ k+1 Mfc+2 Mn 

1 /"ll 

AiS^im-iy. J mk\Jgf e «i~«m<W 71 o^ k+2 n o^n a«r w ... M 

•:/;'-•" (^)r/rr ; , p7U . (9) 



W = J7^-n7Z: L. e ai .-a ro n ai a Mfc+2 n a2 • • • d^n^dx"™ A • • • A (to"" 



- 1)! h h k\ JgT v ~'^ w " w 
where da^...^ is the invariant surface element of Mj and = det(g^). 

From the above discussions, especially the expressions (f|), (||) and (§), we can induce a generalized 
topological current j^ 1 ' " fJ-k which does not vanish only at the zeroes of order parameter field 4>(x), 
and is exactly corresponding to the generalized Winding Number, 

Obviously this tensor current is identically conserved, i.e. 

V w j« = 0, i = l,...,k. 
It is easy to see that j'W"^ are completely antisymmetric tensors. 
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By making use of the 0-mapping theory, we will study the global property of the generalized 
topological current j m " Mfc on the whole manifold G and conclude that j w ' behaves itself like the 
generalized function 5(0). From (JJ) we have 



1 1 1 



which should be looked upon as generalized functions |T8]. Using these expressions the generalized 
topological current (|T0|) can be rewritten as 



J . c c ai ... am 

where C m is a constant 

r — ) ~ A(s m - 1 )(m-2)(m-iy. > m> 2 

and G m (||0||) is a generalized function 

G (\\6\\) = \ Twk^ ' m > 2 
tjmlll0llj \ In 1 10| | , m = 2 * 

Defining general Jacobians J m ' as following 



and by making use of the m-dimensional Laplacian Green function relation|10 



1 47T m/2 

where = ( gjtg^a ) is the m-dimensional Laplacian operator in 0-space, we do obtain the 5- 
function like topological current rigorously 

jUVHk = _5(0)JW-Mfc(^). ( 12 ) 

We find that j^i-Vk ^ o only when = 0, which is just the singularity of j^ 1 "'^ k . In detail, the 
Kernel of the 0-mapping is the singularities of the topological tensor current jW"^ in G. We think 
that this is the essential of the topological tensor current theory and 0-mapping is the key to study 
this theory. 
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To investigate the topological properties of the generalized topological current, we should study 
the total expansion of the 5-function 5(4>). As is well known fl9[| , the 5-function 5(Ni) in curved 
space-time on a submanifold iVj is 

1 



-J n (x - Zi(u , u ))y/g^,d u 



y/g x 

and, by analogy with the procedure of deducing 5(f(x)), since 



(13) 



+00, for <p(x) = 
0, for ${x) ^ 



+00, for x G Ni 
0, for x ^ Ni ' 



(14) 



we can expand the 5-function 5(<j)) as 



5(0) = £Q5(iW 



(15) 



8=1 



where the coefficients c$ must be positive, i.e. q =| q |. From the definition of PVj in (|8|), the Winding 
number can also be rewritten in terms of the parameters v A of Mj as 



_ f Ai—Am a ,n a i...f). n am rl m i) 

Z7T JS, 



Then, by duplicating the above process, we have 



5U)J(-)d m v, 

Si V 



(16) 



Substituting flT5| ) into (0), and considering that only one pi G Sj, we can get 

1 



V 



Si JA^i 



5 n (f - z t (u\ u I ))J[^)^g~d k uJg~d m v. (17) 



where g v = det(^s). Because ^Jg^^fg^d k ud m v is the invariant volume element of the Product 
manifold Mj x Ni, so it can be rewritten as ^fg^d n x. Thus, by calculating the integral and with 
positivity of q, we get 



Ci = zr = ~r~ , 18 

where j3i = \W{\ is a positive integer called the Hopf index [[20f l of 0- mapping on Mj, it means that 
when the point v covers the neighborhood of the zero point pi once, the function covers the 
corresponding region in 0-space $ times, and r/i = signJ(^) Pt = ±1 is the Brouwer degree of </>- 
mapping ]20|]. Substituting this expression of Cj and (|15"D into (|l|), we gain the total expansion of the 
rank-/c topological current 



or in terms of parameters y A = (v 1 , • • •, v m , u 1 ,- ■ •, u k ) 

From the above equation, we conclude that the inner structure of j^ 1 "'^ or j A x'" A k is labelled by the 
total expansion of 5(4>), which includes the topological information j3{ and rji. 

It is obvious that, in (|6|), when u 1 and« 7 (J = 2, k) are taken to be time-like evolution parameter 
and space-like parameters, respectively, the inner structure of j^ 1 '"^ or j A i"' A k just represents I 
(k — l)-dimensional topological defects moving in the n-dimensional Riemann manifold G. The k- 
dimensional singular submanifolds iVj (i = 1, • • -l) are their world sheets. Here we see that the defects 
are generated from where = and, the Hopf indices $ and Brouwer degree rji classify these defects. 
In detail, the Hopf indices /3j characterize the absolute values of the topological quantization and the 
Brouwer degrees rji — +1 correspond to defects while rji — — 1 to antidefects. It must be pointed that 
the relationship between the zero points of the m-dimensional order parameter field <fi and the space 
position of these topological defects is distinct and clear and it is obtained rigorously without tie on 
any concrete model or hypothesis. Furthermore, for the first time we gain the topological charges of 
these defects which are determined by the Winding numbers of the 0-mapping. 



3. The branch processes of the topological defects 

In this section, we will discuss the branch processes of these topological defects. In order to 
simplify our study, we select the parameter u l as the time-like evolution parameter t, and let the 
space-like parameters u 1 = a 1 (I = 2, k) be fixed. In this case, the Jacobian matrices J Al "' Ak (y) 
are reduced to 

y y y y v 

A,B = l,...,(m + l), I j = m + 2,...,n, (20) 

for y A = v A (A < m), y m+1 = t, y m+I = a 1 (I > 2). In the above section, we have studied the 
topological property of the topological defects in the case that the vector order parameter only 
consists of regular points, i.e. ([?[) is hold true. However, when this condition fails, the above results 
will change in some way. It often happens when the zeros of include some branch points, which 
lead to the branch process of topological current. The branch points are determined by the m + 1 
equations 

(f) a (v\- = 0, a = l,...,m (21) 
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and 

4> m+1 (v\--.,v m ,t,a) = J(%) = (22) 

for the fixed a. and they are denoted as (t*,pi). In 0-mapping theory usually there are two kinds of 
branch points, namely the limit points and bifurcation points [2T|, satisfying 



J\pk*, pi )^0 (23) 

and 

J 1 fy\(**i) = 0, (24) 

respectively. In the following, we assume that the branch points (t*,pi) of 0-mapping have been 
found. 

A. Branch process at the limit point 

In order to use the theorem of implicit function to study the branch process of topological defects 
at the limit point, we use the Jacobian J 1 ^) instead of J(^) to discuss the problem. In fact, this 
means that we have replaced the parameter t by v 1 . For clarity we rewrite the problem as 

<f) a (t,v 2 ,---,v m ,v l ,a) = 0, a = l,...,m. (25) 

Then, taking account of the condition(^) and using the implicit function theorem, we have an unique 
solution of the equations (|25|) in the neighborhood of the limit point (t*,Pi) 

t = t{v\a), = v\v 1 ,a), z = 2,3, ...,m (26) 

with t* = t(pj, a). In order to show the behavior of the defects at the limit points, we will investigate 
the Taylor expansion of (p6|) in the neighborhood of (t*,Pi). In the present case, from (|23|) and (|22|), 

we get 

~jjT\(t*,pi) - JJJj \(t*, Pi ) - oo, {M) 

i.e. 

dt l -n 

Then we have the Taylor expansion of ( fffi ) at the point (t*,Pi) 

dttt i / 1 l \ ^ ^ i / i i \ 2 



x -p { 
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± ^~ ^ ^ i /I 1 \ 2 

Therefore 

*-** = ^(^jl(*-*)( vl -P» a ( 28 ) 
which is a parabola in the v 1 — t plane. From (p8|) , we can obtain the two solutions v^(t,a) and 



V( 2 j(t,a), which give the branch solutions of the system (|2l|) at the limit point. If r^rp \ {t*.p t ) > 0, 
we have the branch solutions for t > t*, otherwise, we have the branch solutions for t < t* . The 
former is related to the origin of the topological defects at the limit points. Since the topological 
current of the topological defects is identically conserved, the topological quantum numbers of these 
two generated defects must be opposite at the limit point, i.e. Pitji + — 0. In fact, these two 
cases are just related to the generation and annihilation of defect-antidefect pair. The result (f27j) 
agrees with that obtained by Bray p2| who had a scaling argument associated with point defects 
final annihilation which leads to large velocity tail. 

B. Branch process at the bifurcation point 

In the following, let us consider the case (|24j), in which the restrictions of the system (|2lD at the 
bifurcation point (t*,Pi) are 

J(^)\(t*, Pi ) = 0, J(^)| ( *.«) = 0. (29) 

These two restrictive conditions will lead to an important fact that the dependency relationship 
between t and v 1 is not unique in the neighborhood of the bifurcation point (t*,pi). In fact, we have 

dv \ - JH '\ m) 

which under the restraint ( p9D directly shows that the tangential direction of the integral curve of 
equation ( p0|) is indefinite at the point (t*,Pi). Hence, ( p0|) does not satisfy the conditions of the 
existence and uniqueness theorem of the solution of a differential equation. This is why the very 
point (£*, Zi) is called the bifurcation point of the system (|2T|). 

In the following, we will find a simple way to search for the different directions of all branch curves 
at the bifurcation point. As assumed that the bifurcation point (t*,Pi) has been found from ( |2"T| ) and 
(|2"2"D, the following calculations are all conducted at the value (t*,Pi). As we have mentioned above, 
at the bifurcation point (t*,Pi), the rank of the Jacobian matrix [|^] is smaller than m. In order to 
derive the calculating method, we consider the rank of the Jacobian matrix [|^] is m — 1. The case of 
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a more smaller rank will be discussed in next subsection. Suppose that one of the (m — 1) x (m — 1) 
submatrix Ji(^) of the Jacobian matrix [|^] is 



■hi 



( <P\ 

ft 



V 



im-1 Jjn— 1 



im- 1 
m 



(31) 



and its determinant det Ji(^) does not vanish at the point (t*,Pi) (otherwise, we have to rearrange 
the equations of (]2ID), where 4> a A stands for (dcjf /dv A ) (a — 1, m — 1; A — 2, m). By means of 
the implicit function theorem we obtain one and only one functional relationship in the neighborhood 
of the bifurcation point (t*,Pi) 



v A = f A (v\t } a 



A = 2,3,... ,n 



(32) 



with the partial derivatives 



dv 1 ' Jt 



dv A 
~dt 



A = 2,3, ...,n. 



Then, for a = 1, m — lwe have 



La ±af„,l r2/,l 



(v\r(v\t,a),...,f m (v\t,a),t,a) = 



which gives 



A=2 



dv 1 



A=2 



2^ dyA Jt m 



a = 1, m — 1 



a = 1, m — 1 



(33) 



(34) 



from which we can calculate the first order derivatives of f A : f A and f A . Denoting the second order 
partial derivatives as 

A _ d 2 v A A _ d 2 v A A _ d 2 v A 
fu ~ {dv 1 ) 2 ' ht ~ dv l dt' ftt ~ ~W 

and differentiating fl3"3]) with respect to v 1 and t respectively, we get 



E r A fii = - E Mi/f + E WABf? )fi) - <Pn, a = 1, 2, m - 1 



(35) 



A=2 



A=2 



B=2 



A=2 A=2 



E 



E [4>%fi + Fjaft + E (^ABff)fi] - <f>it, a = 1, 2, m - 1. (36) 



B=2 
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And the differentiation of fl34|) with respect to t gives 

m m m 

E tftf? = - E [20V 4 A + E (<n*/f )f t A ] - 

A=2 A=2 B=2 



1,2, 



m 



(37) 



where 



<f>l 



d 



,2 J>a 



<9 



2lo 



The differentiation of fl34j) with respect to u 1 gives the same expression as fl36|). If we use the Gaussian 



elimination method to the three vectors at the right hands of the formulas (|35|), (|36D and fl37|) , we 
can obtain the three partial derivatives f^, fg. and f£. Notice that the three equations fl35|) , ([36|) 
and (|37j) have the same coefficient matrix Ji(-), which are assumed to be nonzero, and we should 
substitute the values of the partial derivatives and f t f, which have been calculated out in the 
former, into the right hands of the three equations. 

The above discussions do not matter to the last component ^(v 1 , ■ ■ -, v m , t, a). In order to find 
the different values of dv 1 /^ at the bifurcation point, let us investigate the Taylor expansion of 
^(v 1 , ■ ■ -,v m ,t,a) in the neighborhood of (t*,pi). Substituting the existing, but unknown, depen- 
dency relationship ( |3"2| ) into (f) m (v 1 , ■ ■ •, v m , t, a), we get the function of two variables v 1 and t 



F(t,v\a) 



m {v\f{v\t,a),...,r{v\t,a),t,a) 



(38) 



which according to (^) must vanish at the bifurcation point 



F(t*, Pi )=0. 



(39) 



From (|38"D, we can calculate the first order partial derivatives of F(t,v l ,a) with respect to v 1 and t 
respectively at the bifurcation point (t*,pi) 



ar m ap m 

^i = 0r+E <n ft = c + E wtf- 



A=2 



Using (|33|) and (H), the first equation of (p9|) is expressed by 



J (-)\{t*, Pi ) 



- E 

m 

- E 

A=2 



- E - 

A=2 



jm-1 



<,m.— 1 



m 



(40) 
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which, by Cramer's rule, (j3l|) and (0), can be rewritten as 








A=2 



m 



(**,Pi) 



Since 



— detJi(-)| (t . )Pi ) = 0. 



detJiQl^^O 



which is our assumption, the above equation leads to 

dF 



dv 1 



With the same reasons, we can prove that 



dF 
~dt 



\(t*, Pi ) = o. 



\(t*,Pi) = °- 



(41) 



(42) 



The second order partial derivatives of the function F(t, v l , a) are easily to find out to be 

8 2 F 



(dv 1 ] 



+ E + m& + E {^fi)fx\ 



A=2 



B=2 



d 2 F 
dv 1 dt 

d 2 F 



<% + E Wktf + C^/i 4 + mi + E (« fl ) fi] 



A=2 



B=2 



dt 2 



<% + e im t ft A + 4%f£ + E m B ft B )ft 



A=2 



B=2 



which at (t*,pi) are denoted by 

d 2 F 



A 



Then, by virtue of (|3TS|), fl4"T|), ( |4"2"| ) and (T43|), the Taylor expansion of F(t,f 1 ,a) in the neighborhood 
of the bifurcation point (t*,Pi) can be expressed as 



B 



d 2 F 



d 2 F 

\(t*,Pi)i C = -577* \{t*,Pi) ■ 



(43) 



F(t, v\ a) = l ~A{J - p\f + B(v l - p]){t - O + \C{1 - ef 



(44) 



which is the expression of 



Lml „,1 



v m ,t,a) in the neighborhood of (t*,pi). The expression ([Hj) 



shows that at the bifurcation point (t*,pi) 

Aiv 1 - p]) 2 + Wiv 1 - p]){t - f) + c(t - ef 



0. 



(45) 
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Dividing fl45|) by {y 1 — p\) 2 or (t — t*) 2 , and taking the limit t — > t* as well as — > p\ respectively, 



we get two equations 

C(^) 2 + 2B^- 1+ A = 0. (46) 

and 

^lir) + 2B ir + c = °- (47 > 

So we get the different directions of the branch curves at the bifurcation point from the solutions of 
© or §3). 

In order to determine the branches directions of the remainder variables, we will use the relations 
simply 

dv A = ffdv 1 + ffdt, A = 2, 3, n 
where the partial derivative coefficients and / t A have given in (|33|) and (p4|) . Then, respectively 

— = f A + f A — 

dv 1 1 * dv 1 

or 

where partial derivative coefficients ff and f A are given by (|33D and ([34]) . From this relations we 
find that the values of dv A /dt at the bifurcation point (t*,Zi) are also possibly different because 
(J?) may give different values of dv 1 /dt. The above solutions reveal the evolution of the topological 



defects. Besides the encountering of the defects, i.e. two defects encounter and then depart at the 
bifurcation point along different branch curves, it also includes splitting and merging of defects. 
When a multicharged defects moves through the bifurcation point, it may split into several defects 
along these two different branch curves. On the contrary, several defects can merge into one defect 
at the bifurcation point. The identical conservation of the topological charge shows the sum of the 
topological charge of final defects must be equal to that of the initial defects at the bifurcation point, 
i.e., 

X h : 'U. (49) 

/ i 

for fixed j. Furthermore, from above studies, we see that the generation, annihilation and bifurcation 
of defects are not gradual changes, but start at a critical value of arguments, i.e. a sudden change. 

C. Branch process at a higher degenerated point 
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In the following, let us discuss the branch process at a higher degenerated point. In the above 
subsection, we have analysised the case that the rank of the Jacobian matrix [d<p/dv] of the equation 
fl22|) is m — 1. In this section, we consider the case that the rank of the Jacobian matrix is m — 2 
(for the case that the rank of the matrix [d(f>/dv] is lower than m — 2, the discussion is in the same 
way). Let the (m — 2) x (m — 2) submatrix J 2 (-) of the Jacobian matrix [d<f>/dv\ be 



j 2 n 



/ <f>l 



4 \ 



V 



km— 2 
m 



and suppose that det J 2 (j)|(t* )Pi ) 7^ 0. With the same reasons of obtaining (j32D , we can have the 
function relations 



A = 3,4, ...,m. 



(50) 



For the partial derivatives ff, / 2 and f t , we can easily derive the system similar to the equations 
(|33"D and (|34"D, in which the three terms at the right hand of can be figured out at the same time. 
In order to determine the 2-order partial derivatives f^, f^, /^|, f£t an d /if, we can use the 
equations similar to Q55D, ( P5| ) and Q3"7|). Substituting the relations (^Dj) into the last two equations 
of the system (CT), we have the following two equations with respect to the arguments v 1 , v 2 , t,a 



F 1 (v\v 2 ,t,a) = <p m -\v\v 2 ,f 3 (v\v 2 ,t,a),---,f m (v\v 2 ,t,a),t,a) = 



(51) 



\ F 2 (v\ v 2 , t, a) = 4> m (v\ v 2 , / V, « 2 , *, • • • , f m (v\v 2 , t, a),t, a) = 0. 

Calculating the partial derivatives of the function F\ and F2 with respect to v 1 , v 2 and t, taking 
notice of fl50"D and using six similar expressions to ( |iTD and (|42"|), i.e. 



dF, 
dv 1 



\(t*,PiY- 



0. 



dF) 

dv 2 



\(t*,Pi) z 



0. 



d_F) 
dt 



\(t*,Pi) z 



0. 



J = 1,2, 



(52) 



we have the following forms of Taylor expressions of Fi and F 2 in the neighborhood of (t*,Pi) 
F 3 (v\v 2 , t, a) « A^ 1 - p, 1 ) 2 + ^(^ - p])(v 2 - p 2 ) + A^v 1 - p}) 
(t _ t *) + A 3i (v 2 - p 2 ) 2 + A j5 (v 2 - p 2 )(t - t*) + A j6 (t - ff = 

J = 1,2. 



(53) 



In the case of Aj\ 7^ 0, A,- 4 7^ 0, by dividing (^) by (£ — t*) 2 and taking the limit t — > t*, we obtain 



two quadratic equations of and 



j ~dt 



A j2 



dv 1 dv 2 
dt dt 



* dvl , 



+ A j5 - + A je 



(54) 
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Eliminating the variable dv 1 /dt, we obtain a equation of dv 2 /dt in the form of a determinant 



A u A 12 Q + A 23 A U Q 2 + A 15 Q + A 16 

A u A 12 Q + A 13 A lA Q 2 + A 15 Q + A 16 

A 2 i A 22 Q + A 23 A M Q 2 + A 25 Q + A 26 

A 21 A 22 Q + A 23 A 24 Q 2 + A 25 Q + A 



26 



= 



(55) 



where Q = dv 2 /dt, which is a Ath order equation of dv 2 /dt 

,dv 2 .4 ,dv 2 ,n ,dv 2 ,o ,dv 2 , 
«o(-r-) 4 + ) 3 + a 2 (^-) 2 + a 3 (—) + a 4 = 0. 

at at at at 



(56) 



Therefore we get different directions at the bifurcation point corresponding to different branch curves. 
The number of different branch curves is four at most. If the degree of degeneracy of the matrix [|^] is 
more higher, i.e. the rank of the matrix [|^] is more lower than the present (m — 2) case, the procedure 
of deduction will be more complicate. In general supposing the rank of the matrix [|^] be (m — s), 
the number of the possible different directions of the branch curves is 2 s at most. Comparing with 
the above subsection, the solutions in the subsection also reveal encountering, splitting and merging 
of the defects along more directions. 

At the end of this section, we conclude that there exist crucial cases of branch processes in our 
topological defect theory. This means that a topological defect, at the bifurcation point, may split 
into several topological defects along 2 s different branch curves with different charges. Since the 
topological current is a conserved current, the total quantum number of the spliting topological 
defects must precisely equal to the topological charge of the original defect i.e. 

2 s 

for fixed i, P^r)^ stands for the total topological charge of defects along the j'-th branch curve at the 
bifurcation point p^. This can be looked upon as the topological reason of the defect splitting. Here 
we should point out that such splitting is a stochastic process, the sole restriction of this process is 
just the conservation of the topological charge of the topological defects during this splitting process. 
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